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ABSTRACT 

We present a backward approach for the interpretation of the evolution of the near-infrared 
and the far-infrared luminosity functions across the redshift range 0 < 2 < 3. In our method, 
late-type galaxies are treated by means of a parametric phenomenological method based on 
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PEP/HerMES data up to 2 ~4, whereas spheroids are described by means of a physically moti¬ 
vated backward model. The spectral evolution of spheroids is modelled by means of a single-mass 
model, associated to a present-day elliptical with K —band luminosity comparable to the break of 
the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr 
across the redshift range 1 < z < 5. The key parameter is represented by the redshift 20.5 at 
which half proto-spheroids are already formed. A statistical study indicates for this parameter 
values between 20.5 = 1.5 and 20.5 = 3. We assume as fiducial value 20.5 ~ 2, and show that this 
assumption allows us to describe accourately the redshift distributions and the source counts. By 
assuming z 0 .5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources 
observed at 2 > 2 can be explained as progenitors of local spheroids caught during their forma¬ 
tion. We also test the effects of mass downsizing by dividing the spheroids into three populations 
of different present-day stellar masses. The results obtained in this case confirm the validity of 
our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model 
which describes the evolution of galaxies at the break of the present-day early type A'-band LF. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: luminosity function, mass 
function — infrared: galaxies 


1. Introduction 


Achieving a complete understanding of the ori¬ 
gin of the local dichotomy of spheroids and disc 
galaxies has been one of the main objectives in 
astrophysics during the past years, as well as ob¬ 
taining an accurate measure of the star formation 
history over cosmic time. 

From the pioneering IRAS satellite results in 
the local Universe (e.g. Soifer et al. 1987) and 


the detection of a cosmic infrared background as 
energetic as the optical/near-IR background (e.g. 
Puget et al.| 1996), it is now well established that 


most energy radiated by newly formed stars is 
heavily absorbed by dust and re-emitted in the 
infrared (IR) band. In the last decade, the ISO 
and the Spitzer satellite individually detected IR 
sources up to 2~1 (e.g. 


pioni et al. 

2002 

) and 2 

2004 

Shupe et al. 2008) 


Elbaz et al. 

1999 

Grup- 

~2 (e.g. 

Papovich et al. 


their capabilities at far-IR wavelengths (i.e. where 
the dust reprocessed emission peaks) were still 
strongly limited by source confusion. 

Nowadays, the Herschel Space Observatory 
(Pilb ratt et al.|201 0) . las allowed to properly mea¬ 
sure the IR luminosity function of galaxies up to 
2~4 ( Gruppioni et al.||2013 hereafter GPR13], see 
also Magnelli et al. 2013), thanks to its mirror of 


3.5-m and an observing spectral range between 60 
and 670 pm. The derived IR luminosity density 
( pm ) confirms the Spitzer 24-pm based results 
(e.g. Caputi et al. 2007, Magnelli et al. 2009 


Rodighicro et al.||2010) up to 2 ~2, revealing that 


the IR luminosity density increases steeply from 
2 = 0 up to 2 ~ 1 , then flattens between 2 ~ 1 and 
2^3, to decrease at 2 >3. 


In the recent paper by Burgarella et al. (20131, 
the IR Herschel luminosity function derived by 
GPR13 has been combined with the luminosity 
function in the far-UV from |Cucciati et a l. (20121, 
in order to achieve an estimate of the redshift evo¬ 
lution of the total (far-UV + IR) star formation 
rate density (SFRD). The SFRD is always domi¬ 
nated by the IR emission, whereas the UV contri¬ 
bution increases steeply from 2 = 0 up to 2^2.5, 
where it flattens and settles on a plateau up to the 
highest redshift sampled by the survey (2 ~ 3.6). 
This suggests that the peak of the dust attenua¬ 
tion, occurring at 2 ~ 1, is delayed with respect to 
the SFRD plateau (2 ~ 2 — 3), derived from the 
far-UV. 


The accurate determination of the star-formation 
history up to 2~4 has rendered particularly ur¬ 
gent the issue of explaining theoretically this be¬ 
haviour, as well as understanding how the ob¬ 
served SFRD evolution is linked to the galaxy 
formation process. 


Within the current ‘concordance’ cosmologi¬ 
cal paradigm, which sees a A-Cold dark Mat¬ 
ter (ACDM)-dominated universe, the formation of 
structures is hierarchical, since small dark matter 
(DM) halos are the first to collapse, then inter¬ 
act and merge to assemble into larger halos (e.g. 
Lacey et al.||2008 Fontanot et al.|[2009 ). 


In cosmological ACDM-based semi-analytical 
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models (SAMs) for galaxy formation, the most un¬ 
certain assumptions concern the behaviour of the 
baryonic matter. In the first classical SAMs im¬ 
plementations, baryonic matter was assumed to 
follow the DM in all the interaction and merg¬ 
ing processes and spheroid galaxies were formed 
from several merging episodes of smaller sub-units, 
with the most massive galaxies as the last systems 
to assemble. Indeed, more recent results suggest 
that spheroids may form over a wide redshift range 
by means of both galaxy mergers and disk insta¬ 


bilities (e.g. De Lucia et al. 2011 Porter et al. 


2014). In particular, the latter process should be 


of fundamental importance if the vast majority of 
z~2 quasi-stellar objects is associated with unper¬ 
turbed systems, as indicated by recent results (e.g. 
Schawinski et al. 20121. 

Although accounting for the recent updatings, 
SAM galaxy formation models still encounter 
difficulties to reconcile properly with the well- 
established anti-hierachical evolution of spheroids, 
a feature which is known as “downsizing” and 
which is supported by a large amount of obser¬ 
vational data ( Matteucci|| 1994 Cowie et al.||1996 


Kodama et al.|2004[ but see also |Marchesini et al.| 
20091). 


Moreover, ACDM models highlight difficulties 
for reproducing the evolution of IR galaxies (in 
particular the bright-end IR luminosity function 
at high-z, see e.g. |F ontanot fe Somerville || 20TH 
Niemi et al., |20121 Gruppioni, Calura, Pozzi et al. 
in prep.). The basic problem of the dominance of 
a merger-driven evolution is the duration of the 
star-format ion in spheroids, which extends over 
times much longer than those indicated by the ba¬ 
sic stellar population diagnostics of local ellipticals 
( Matteucci||1994 Renzinif~2006 ). 

Doubts on the dominant role of mergers in driv¬ 
ing the star formation histories of galaxies is also 
supported by the observation of systems with reg¬ 
ular kinematics and showing little signs of major 


mergers ( 

Cresci et al. 2009 

Forster Schreiber et al. 

2009 

Gnerucci et al. 

2011 

), as well as the tight 


scaling relations of star-forming galaxies, such as 


the mass-SFR plane ( 

Rodighiero et al. 

2011 

and 

the fundamental metallicity relation 

Mannucci 


et al. 2010 Lara-Lopez et al. 2010), apparently 


incompatible with the stochastic nature of galaxy 
mergers (see e.g. Genzel et al. 20081. 


In Dekel et al. (2009), a new scenario is pro¬ 


posed where galaxy formation is not mainly 
merger-driven, but occurs instead via a ‘cold 
stream accretion’ process, which represents the 
dominant mode for galaxy formation. Such a sce¬ 
nario is able to account for several galaxies prop¬ 
erties at high-z (i.e. gas-richness and extended 
rotationally supported disks). In the proposed 
model, the streams responsible for the creation 
of rotationally-supported discs could also sup¬ 
ply the turbulence necessary for the formation 
of spheroids. 



Cook et al. 

(2009 

), building upon the model by 

Granato et al. 

(2004 

), proposed an alternative see- 


nario for the origin of the local dichotomy between 
spheroids and discs in galaxies, in the context of 
the hierarchical build up of dark matter (DM) 
haloes. Noting that the growth of DM haloes oc¬ 
curs in two phases (e.g. Zhao et aLf 2003), namely 
a fast collapse one, featuring a few major merger 
events, followed by a late, quiescent accretion onto 
the halo outskirts, they proposed that the first 
and second phases of halo growth drive two dis¬ 
tinct modes for the evolution of baryonic matter, 
favoring the development of the spheroidal and 
disc components of galaxies, respectively. As a re¬ 
sult, in that semi-analytic model, the spheroidal 
component develops rapidly and at early times 
(z 2), mimicking the classic ‘monolithic model’ 
(e.g. Larsonp~976 Matteucci||1994 ) 


Although a similar picture can account for 
the duality between spheroids and discs, obtain¬ 
ing a full, theoretical understanding of the basic 
properties of galaxies as traced by current multi¬ 
wavelength surveys over cosmic time, is still a 
challenging task. 

At present, the lack of a complete theoreti¬ 
cal picture for galaxy formation encouraged sev¬ 
eral attempts to reproduce the basic observational 
properties of galaxies over cosmic time via ‘hy¬ 
brid’ models, which generally match simple phys¬ 
ical models to parametric approaches. A few at- 


tempts (i.e. 

Silva et al. 

2005 

Lapi et al. 

2011 

Cai et al. 

2013 

Bonato 

et al. 

2014 

) have been 


able to provide a reasonable description of the 
multi-wavelength properties of galaxies, from the 
near-IR up to sub-mrn wavelengths. However, the 
treatment of the chemical evolution of dust grains 
in such systems has often been simplistic, i.e. as¬ 
suming a constant carbon-to-silicate ratio or dust- 
to-gas ratios simply scaling with metallicity. 
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In the work of Schurer et al. (2009) detailed 


predictions on the evolution of the dust chemi¬ 
cal composition in galaxies of different morpho¬ 
logical types have been matched to the spectro- 


photometric code GRASIL (Silva et al. 19981, 


which takes into account dust reprocessing, to 
analyse the evolution of the resulting spectral en¬ 


ergy distributions (SEDs). Schurer et al. (2009) 


showed that a detailed treatment of the dust chem¬ 
ical composition is particularly important when 
the proto-spheroids SEDs are taken into account, 
since dust assumptions tuned to match the dust 
depletion pattern of the Milky Way turned out as 
particularly inadequate to describe their observ¬ 
ables. 

In this paper, we present a new backward 
phenomenological approach where the spectro- 
photometric evolution of dust in proto-spheroids 
is calculated self-consistently, i.e. from idealised 
star formation histories able to account for a large 
set of observational data at various redshifts, in¬ 
cluding their dust content and local abundances 


ratios (Calura et al. 2008 [hereafter CPM08]; 
Pipino et al. 20111 and the evolution of their SEDs 


(Schurer et al. 2009). Beside our physically mo¬ 


tivated ‘backwards’ model for the evolution of 
proto-spheroids, the evolution of late-types galax¬ 
ies is also taken into account through a parametric, 
phenomenological approach. 

Our aim is to reproduce a variety of multi¬ 
wavelength observables, including source counts, 
redshift distributions and luminosity functions at 
various redshifts and from the near-IR to the far- 
IR. 


The basic questions we aim to address are to 
which extent the star-forming galaxies detected 
by Herschel at high-z, could be considered as 
the progenitors of local early galaxies (i.e. proto¬ 
spheroids), and which quantitative constraints can 
be achieved from current observations regarding 
the dominant epoch for spheroid formation. 

The paper is structured as follows: in Section 

2 we describe how the model deals with the late 
type and proto-spheroids populations, in Section 

3 we compare the far-IR and iGband observables 
with the model predictions, and in Section 4 we 
discuss the formation redshift and mass buildup 
of the proto-spheroids population. 

All magnitudes reported in this paper are in the 


AB system (Oke & Gunn 1983). Throughout this 
paper, we adopt a flat A CDM cosmology with 
H 0 = 71 kms -1 , = 0.27 e D A = 0.73. 


2. Model description 


The basic idea behind our approach is similar to 
the one of previous ‘hybrid’ models used to inter- 


pret far-IR and optical data (see 

Silva et al. 2005 

Lapi et al. 

2011 

Cai et al. 2013 

1, where galaxies 


are divided into two main classes, i.e. early-type 
and late-type galaxies, and each class treated in a 
different way. 


The evolution of late-type galaxies is followed 
by means of a parametric, classical ‘backward’ ap¬ 
proach, which consists in evolving a local lumi¬ 


nosity function (LLF) back with redshift (Gruppi- 


oni et al. 2011). The observables are reproduced 


by evolving the LLF of different populations, each 
characterized by a set of Spectral Energy Distribu¬ 
tions (SEDs), in density and luminosity. This ap¬ 
proach has been traditionally used to interpret far- 


IR source counts (e.g. 

Pozzi et al. 

2004 

Valiante 

et al. 

2009, |Franceschini et al. 2010 ( 

druppioni 

et al. 

2011 

Bethermin et al. 2011). 


The evolution of proto-spheroids is followed by 
means of a model in which the chemical com¬ 
position of dust is calculated in detail and self- 
consistently on the basis of its star formation his¬ 
tory. In this model, the evolution of the SED of the 
proto-spheroid is followed as a function of time. 


The main reason for developing this approach 
derives from the need of finding a model able to 
account for data in various bands. Up to now, 
infrared backward evolutionary models, although 
accurately accounting for the dust far-IR emis¬ 
sion (linked to the ongoing star-formation, hence 
to young stellar populations), have neglected the 
need to account for other observables related 
to evolved stellar populations, such as the K- 
band data, which trace the oldest stellar popu¬ 
lations. Data from IL-band up to far-IR wave¬ 
lengths can be simultaneously reproduced only 
if a self-consistent model that accounts for the 
time-dependent galaxy SED is considered. 

The description of our method can be sub¬ 
divided into three steps. First, we divide far- 
IR sources into late-type galaxies and ‘potential’ 
proto-spheroids on the basis of their SED and, 
more importantly, on the basis of the evolution 
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Fig. 1.— Evolution of the comoving number 
density of the PEP sources up to 2 ~4 (taken 
from Fig. 18 of GPR13). Black filled dots 
with error-bars within the ±lcr uncertainty re¬ 
gion (grey filled area) represent the total sample. 
The coloured shaded areas represent the contribu¬ 
tion from different populations, highlighted as fol¬ 
lows: green: spiral; dark-green: SF-AGN(Spiral); 
cyan: starburst; orange: SF-AGN(SB); magenta: 
AGN2; blue: AGN1. The dark-green dashed 
line corresponds to the sum of the spiral+SF- 
AGN(Spiral) populations; the red dashed line 
corresponds to the sum of the starburst+SF- 
AGN(SB)+AGN2+AGN1 populations. 


(positive or negative) of their comoving number 


density as a function of redshift z (see Sect. 2.1|. 


This approach is different from that used in other 


works, such as Cai et al. (2013), where IR galax¬ 


ies are classified according to their redshift, with 
high- 2 : (z>1.5) star-forming sources generally in¬ 
terpreted as proto-spheroids, and sources at lower- 
2 ( 2 <1.5) as late-type galaxies. Secondly, we nor¬ 
malize the proto-spheroid galaxy number density 
according to the local A-band LF. Finally, we use 
our physical, chemo-spectrophotometric SED evo¬ 
lution (Schure r~et al.|20 09) to calculate the source 
counts and redshift distributions for the proto¬ 
spheroid population. The proto-spheroids forma¬ 
tion redshift is left as the only free parameter; 
for this quantity, we explore an appropriate range 
in order to simultaneously reproduce the A'-band 
and far-IR observables. 

The main difference between the approach used 
in this paper and previous ‘hybrid’ models (Silva 
et al. 2005, Lapi et al. 2011) is that in the 
latter, the abundance of spheroids as a function 
of redshift is calculated by means of a simplified 
treatment of the merging history of dark mat¬ 
ter halos, described by the |Press fe Schechter] 
(1974) formalism. In this paper, the statistics 


of spheroidal galaxies are determined purely on 
a phenomenological basis, i.e. from their present- 
day AT-band LF, with their backwards evolution 
describing their abundance at redshift 2 > 0. 


2.1. IR populations 

Our IR sources division into late-types and 
proto-spheroids is based on the results of GPR13. 
In that paper, the IR luminosity functions (LFs) 
of galaxies and active galactic nuclei (AGN) up to 
z~4 were presented. The IR LFs were derived by 
exploiting the deep and extended data-sets of the 
Herschel GTO PACS Evolutionary Probe (PEP, 
Lutz et al. 2011). Different cosmological fields at 
different depths were included, from the shallow 
COSMOS field (3<r depth of 10.2 mJy at 160 pm) 
up to the pencil-beam GOODS-S field (3cr depth 
of 2.4 mJy at 160 pm). Far-IR Herschel data from 
the PEP survey (at 70, 100 and 160 pm) were ex¬ 
ploited together with data in the sub-mm band 
from the HerMES survey (Oliver et al.| 2010 at 


250, 350, and 500 pm), and taking advantage of 
the extensive multi-band coverage. 

In GPR13, the sources were classified on the ba- 
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sis of an accurate SED-fitting by comparing data 
with semi-empirical models from different libraries 


In detail, spiral galaxies dominate the global 
density of galaxies in the local Universe (followed 


(Polletta et al. 

2007 

Rieke et al. 

2009 

Gruppi- 

by the SF-AGN(Spiral) systems), but their num- 

oni et al. 

2010 

. According to the best-fit tem- 

ber density at 2~1 already drops by nearly one 


plates, sources were divided into five different pop¬ 
ulations: spiral and starburst galaxies (i.e. galax¬ 
ies with no evidence of an active galactic nucleus, 
and characterized by < T dust >~20 and 40-50 
K, respectively), AGN1, AGN2 and star-forming 
(SF)-AGN. The latter class represents the largest 
fraction (~48 %, see also Sajina et al. 2012) of 
the PEP selected sources, and includes systems 
whose SEDs were fitted by means of empirical 
templates of nearby ULIRGs/Seyferts contaninig 
an obscured or low-luminosity AGN. However, the 
obscured/low-luminosity AGN does not dominate 
the energetic budget in these sources; its presence 
is detected mostly in the narrow mid-IR range, 
where the flux of the host galaxy presents a min¬ 
imum. Spiral galaxies constitute the second more 
numerous class (~38 %), followed by starburst 
galaxies representing ~7 % of the far-IR selected 
sources, and lastly by the AGN1 and AGN2 (only 
3 and 4 %, respectively). 

Given the large fraction of sources presenting a 
SF-AGN SED and the large heterogeneity of the 
templates in this category, SF-AGNs in GPR13 
were further divided in two classes, depending 
on their far-IR/near-IR colours (e.g. Si 00 /Si, 6 ): 
specifically, galaxies fitted by cold Syefert2/1.8 
templates were classified as SF-AGN (Spiral) with 
far-IR/near-IR colour typical of Spiral/late type 
sources, while sources showing the presence of a 
warmer dust component, such as NGC 6240, were 
classified as SF-AGN (SB), with colours typical of 
starburst galaxies. 

In Fig. [l] the comoving number density of 
the different populations has been reported, taken 
from Fig. 18 of GPR13 (top panel). The number 
density in GPR13 were computed by integrating 
the modified Schechter functions that best repro¬ 
duced the different populations down to luminos¬ 
ity L ~ 10 s A®. As already noted by GPR13, the 
different populations behaves in significantly dif¬ 
ferent and rather opposite ways: while the comov¬ 
ing number density of spiral and SF-AGN(spiral) 
strongly increases going from the high -2 to the lo¬ 
cal Universe, the number density of starbursts, SF- 
AGN(SB), AGN1 and AGN2 sources significantly 
decreases with decreasing redshift. 


order of magnitude. Starbursts, SF-AGN(SB), 
AGN1 and AGN2, instead, dominate the number 
density at z > 1.5, whereas their contribution to 
the local number density is negligible. We caution 
the reader that the described trends of the number 
densities, in particolar that referred to the spiral 
galaxies, are not in contrast with the observed de¬ 
crease of the star-formation density from 2 ~ 1 to 
2 ~0 (i.e. Hopkins & Beacom 2006), directly pro¬ 


portional to the IR luminosity density (reported 
in Fig. 18 of GPR13, lower panel) and not to the 
number density. 

In this work, we interpret the evolution of the 
IR comoving number density of starbursts, SF- 
AGN(SB), AGN1 and AGN2 sources, progres¬ 
sively decreasing from z ~ 3 to z = 0, as the 
transition of all these systems to a population of 
passive objects. In other words, in our scheme, we 
regard all these classes of objects as the progeni¬ 
tors of local spheroids. In our formalism, the evo¬ 
lution of the SED of this class of objects is taken 
into account by means of the physical model de¬ 
scribed in Sect. 12.21 

On the other hand, on the basis of their SED- 
classification and on the positive evolution of their 
comoving number density as a function of the cos¬ 
mic time, we classify Spiral and SFR-AGN(Spiral) 
sources as late-type galaxies. As already intro¬ 
duced in Sect. [2j for these galaxy classes, we do 
not consider any physical model, but we adopt 
the same parametric LF evolution as reported in 
GPR13. More specifically, to describe the LF evo¬ 
lution of these populations at different cosmologi¬ 
cal times, we consider both the LLF (computed as 
a modified Schechter function) and the evolution 
of the L star and <F stal . parameters as presented in 
GPR13. 

2.2. The proto-spheroid model 

A long-lasting debate is still open on the for¬ 
mation of spheroids. The dominant merger-driven 
scenario, besides remarkably successes, encounters 
several difficulties in reproducing some basic prop¬ 
erties related to the scaling properties of spheroids 


and their evolution ( 

Matteucci 1994 

Renzini 2006 

Ciotti 2009 Pozzetti et al. 2010 

). For this reason, 
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the old idea (i.e. Larson 1976) that large ellip¬ 


tical galaxies formed their stars at high redshift 
in a huge burst of star-formation is nowadays be¬ 
ing revisited (see Dr iver et al.| [2013 and references 
therein). 

The model we adopt for proto-spheroids is 
based on this idea, and is similar to the classic 
‘monolithic model’ proposed by various authors 
(e.g., Larson 1976 Matteucci 1994). More de¬ 


tailed descriptions of the model can be found in 
CPM08 


v = 15 Gyrs^ 1 and a Salpeter initial mass func¬ 
tion in the mass interval 0.1-100 M 0 . Under these 
conditions, the time of the onset of the galactic 
wind is ~ 0.7 Gyrs. 

In the model, the evolution of the chemical com¬ 
position of the ISM is followed in detail and as a 
function of time. The model takes into account 
also how the interstellar dust composition varies 
as a function of time, instead of assuming a non¬ 
evolving dust composition similar to that of the 


( 2011 ). 


Schurer et al. 

(2009 

and 

Pipino et al. 

Milky Way (seejDraine & Li 2007) 


The model described here is designed to re¬ 
produce an elliptical galaxy of present-day stellar 
mass M* ~ 10 11 M 0 , which corresponds to the 
stellar mass at the break of the local early-type 
stellar mass function ( Baldry et al.|2004 |. Its ba¬ 
sic assumptions are the following. 

In our scheme, proto-spheroid galaxies form as 
the result of a rapid collapse of a homogeneous 
sphere of primordial composition, generating an 
intense star-formation event (see also |Pipino et al. 
2002 Pipino et al. 2005). After the initial collapse, 


the galaxy is allowed to evolve as an ‘open box’ 
into the potential well of a dark matter halo. 

The system is assumed to accrete gas via an 
infall episode; the infall rate as function of time t 
can be expressed as: 


(G)inf OC exp-t/Tinf. 


Cl) 


The quantity r,; rj / determines the timescale of the 
collapse; in this paper, we adopt Tj n / = 0.4 Gyrs 
( Pipino et al.||20lT| ). 

The SFR (expressed in M 0 /yr) is calculated as: 


ip(t) = vM gas (t), 


( 2 ) 


i.e. it is assumed to be proportional to the gas 
mass M gas via the star formation efficiency v, ac¬ 


cording to the Schmidt (1959) law. 


The star formation is assumed to halt as the 
energy of the ISM, heated by stellar winds and su¬ 
pernova (SN) explosions, exceeds the binding en¬ 
ergy of the gas. At this stage, a galactic wind 
triggers, sweeping away all the residual gas. From 
this point on, the galaxy is assumed to evolve pas¬ 
sively. 


As described in detail in Schurer et al. (20091, 


Also, our model includes a detailed treatment 
of the rate of Type la SNe, which are assumed 
to originate from white dwarfs in binary systems, 
as well as a computation of their energetic feed¬ 
back to the ISM (Pipino et al.| |2005 ). These SNe 


continue exploding even after star formation has 
been quenched by the galactic wind. For this rea¬ 
son, they may have an important role in keeping 
the galaxy free of gas and passively evolving. In¬ 
deed, these events occur in a medium already hot 
and rarefied, where radiative cooling is inefficient 
and radiative losses are negligible, so that most of 
their blast wave energy can be transferred as ther¬ 
mal energy into the interstellar medium. Actually, 
in the context of our model, no additional source of 
feedback is required to avoid subsequent episodes 
of star formation. It is however worth noticing 
that semi-analytical models, which, at variance 
with respect to our simplified treatment, do in¬ 
clude a computation of the gas accretion from the 
intergalactic medium, and in recent years begun to 


take into account the effects of SNe la (e.g. Yates 


et al. 2013 Gargiulo et al. 2014), still indicate that 


‘radio-mode’ AGN feedback is required in order to 
avoid low-redshift star formation episodes in ellip¬ 
ticals. 

The chemical evolution model for proto-spheroids 
of CPM08 has then been used by [Schurer et al.| 
(2009) to calculate how the theoretical SED is 


expected to vary as a function of time, once the 
chemical evolution of interstellar dust is taken 
into account in a self-consistent way. To this 
purpose, the spectro-photometric code GRASIL 


(Silva et al. 1998) has ben used. The unique fea- 


in our model we assume a star formation efficiency 


ture of GRASIL is that it allows to track the lumi¬ 
nosity evolution of composite stellar populations 
by taking into account the effects of interstellar 
dust. 

In this paper, the spectral evolution of the 
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proto-spheroid population is modelled by means of 
the theoretical SEDs calculated by ISchurer et alTl 
(2009). 

No attempt is made to model composite galax¬ 
ies, i.e. objects including both bulges and disks. 
These systems are particularly complex since they 
may present composite star formation histories 
and geometries for their dust and stellar compo¬ 
nents in principle different than those assumed in 
this paper. 


3. Results 


In this section, we present the model pre¬ 
dictions (in terms of source counts, redshift- 
distribution and luminosity function), together 
with the corresponding observables. As reference 
bands we choose the A'-band, probing the old 
stellar populations, and the far-IR 160 (im band, 
probing the still-forming stellar populations. 


For the far-IR we use as reference the PEP 
survey. In this survey, sources are selected at 
100 /irn and 160 pm (also at 70 /mi, but only in the 
GOODS-S field). We consider the 160 /tm sample 
as our reference sample, to allow a direct compar¬ 
ison with the results of GPR13, also based on the 
160 /im PEP catalogue. Moreover, the 160 /tm 
sample reaches slightly higher redshifts. 


As already mentioned in Sect. [5] we normalize 
the number density of late-type galaxies according 
to the IR LF as determined by GPR13. The LF in 
the AT-band as a function of redshift is calculated 
on the basis of appropriate colour transformations 
from the best-fitting SED templates. 

At x = 0, the assumed LF of spheroids is the lo¬ 
cal luminosity function of early-type galaxies from 


Kochanek et al. (2001); at redshift greater than 


zero and in the IR, the proto-spheroid LF is deter¬ 
mined on the basis of the evolution of their SED, 
calculated self-consistently by means of the chemo- 
spectrophotometric model described in Sect. |2.2| 
No evolution is assumed for the faint-end slope a. 
At redshift greater than z s , corresponding to the 
beginning of the proto-spheroids formation epoch, 
the normalisation is set to zero. 


3.1. Source counts and redshift distribu¬ 
tions 

At any given flux S*, we calculate the differ¬ 
ential number counts (per unit solid angle) of 
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Fig. 2.— Top panel: Merit function as defined 
in Sect. |3.2| as a function of zq . 5. The red dashed, 
green dashed and black solid lines are the merit 
function calculated considering only the A'-band 
and 160 pm redshift distributions, the A'-band and 
160 pm source counts and the four observables al¬ 
together, respectively. Middle panel: 160 pm red¬ 
shift distribution obtained considering zq .5 = 1.2. 
Bottom panel: 160pm source counts obtained con¬ 
sidering Z 0.5 = 3.5. 
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Fig. 3.— Upper panel: 160 p m differential extra- 
galactic source counts normalized to the Euclidean 
slope. Herschel data from PEP (Berta et al. 2011) 
are represented by red filled circles (GOODS-N), 
blue filled squares (Lockman Hole), green trian¬ 
gles (COSMOS) and cyan pentagons (GOODS-S). 
The grey shaded area represents the 1 — a uncer¬ 
tainty region of the Spitzer data (|Dole et al. 2004; 
Frayer et al. 2006; Bethermin et al. 2010). The 


Herschel data from the combined PEP/GOODS- 


H map (Magnelli et al. 2013) are represented by 


yellow dots (GOODS-S ultradeep) and yellow pen¬ 
tagons (GOODS-N/S deep). Lower panel: Red- 
shift distribution of the PEP sources with 160 /irn 
flux >20 mJy (grey histogram, bin=0.2). The 
model curves for the different populations are plot¬ 
ted in different colours (black solid: total; dotted 
green: spiral; dot-dot-dot-dashed dark-green: SF- 
AGN(Spiral), long-dashed dark-green: spiral+SF- 
AGN(Spiral); dashed red: proto-spheroids). For 
the model results, a piecewise constant continuous 
formation for proto-spheroids has been assumed as 
described in Sect. 13221 


Fig. 4.— Upper panel: Observed and pre¬ 
dicted A'-band differential extragalactic source 
counts. Observational data are from the 
Durham public compilation (website: http://star- 
www.dur.ac.uk/~nm/pubhtml/counts), from 
Madau & Pozzetti (2000) and from the Ul- 


traVISTA survey (95% complete to AT=23.8, 


McCracken et al. 2012), and are shown as yel¬ 


low circles, cyan triangles and pink triangles, 
respectively. The curves represent the same 
model predictions as in Fig. [3] Lower panel: 
Redshift distribution of the A'-band selected 
sources from two surveys. Top: K20 survey 


(Mignoli et al. 2005) at a magnitude limit of 


A'~21.9 (K[Vega] = 20) and over an area of 
55 arcmin 2 (bin=0.05). The curves represent the 
model predictions shown in Fig. [3] 
















































proto-spheroids galaxies from the integral: 




Fig. 5.— Theoretical differential proto-spheroids 
source counts in various redshift bins. Upper 
panel: A'-band counts. Lower panel: 160 /jm 
counts (normalised to the Euclidean slope). A 
piecewise constant continuous formation for proto¬ 
spheroids has been assumed as described in 
Sect. 13.21 


dN 

dSdn 


r^Arr^c* \ , dL(S*, z) dV 
(S*)=/ *[ £(S *, 2 ),z]^-—* 


( 3 ) 

where <F(A, z) is the luminosity function calcu¬ 
lated at the redshift z, dV/dzdVl is the comov¬ 
ing volume per unit of solid angle and z max is the 
maximum redshift where the proto-spheroids LF 
is sampled (i.e. corresponding to the formation 
redshift.) The relation between the rest-franre lu¬ 
minosity L and the observed flux S* at a given 
redshift 2 is: 


logL(S *, z) = logS*- 1-2 log DL+const+E(z)+K(z) 

(4) 

where D^ is the luminosity distance, E(z) is the 
evolutionary correction, K{z) is the K-correction 
and const a constant term (Poggianti ( |1997 |; Shi- 
masaku & Fukugita (1998])). 


Similarly, we calculate the source redshift dis¬ 
tribution (i.e. number of sources above a given 
flux limit Su m , in a redshift range Z\ < z < z-i 
and per uniti of solid angle) from the integral: 


N(> Slim) = 




<t>[L(S,z),z\dL(S,z)—dz 

( 5 ) 


where L m in(Si im ,z) is the rest-frame minimum lu¬ 
minosity corresponding to the flux limit Su m at a 
given redshift z. 


3.2. A continuous formation for proto¬ 
spheroids 


In this work, we assume a continuous forma- 
tioiQof proto-spheroids, described by the function 
'(A and occurring throughout an extended redshift 
range, which starts at redshift z s and ending at z e . 


The local number density of spheroids is de¬ 
termined by the local AT-band LF of |Kochanek 


et al. (2001). For the sake of simplicity, we as¬ 
sume that ^ is a mass-independent function over 


1 Throughout this work, the proto-spheroids in formation are 
those which have started forming stars at a given redshift. 
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Fig. 6. - Observed and predicted AT-band rest-frame luminosity function in eight different redsliift bins. The 
model curves for the different populations are plotted in different colours (black solid: total; dotted green: 
spiral; dot-dot-dot-dashed dark-green: SF-AGN(Spiral), long-dashed dark-green: spiral-f-SF-AGN(Spiral); 
dashed red: proto-spheroids). Observational data in the local Universe are from Kochanek et al. 2001 (grey 
squares). In the local Universe plot, the red and green shaded areas represent the 1-cr uncertainty regions of 
the observed LF for early-type and late-type galaxies, respectively. Data points at higher redshift are from 
Pozzetti et al.|[2003 (open squares) and from Cirasuolo et al.pOlO (grey circles). 
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a given redshift range Az = z s — z e . By means 
of this function, the global spheroid population is 
divided in various sub-populations, each one born 
at a different epoch, thus characterised by a differ¬ 
ent age at the present-day. Each sub-population 
is evolved backwards in luminosity by means of a 
single spectro-photometric model calculated for a 
galaxy of present-day stellar mass of ~ lO n M 0 , 
corresponding roughly to the break of the present- 
day early type A'-band LF. 

We believe that z e = 1 can be regarded as a rea¬ 
sonable value for the end of the formation epoch 
of M> 10 11 Mq spheroids. This value is consistent 
with their stellar populations diagnostics, glob¬ 
ally indicating that the bulk of the most-massive 


early types were already in place at z = 1 (Ren- 


observables represented by the source counts and 
redshift distribution (at the flux limit of the COS¬ 
MOS field) at 160 //m and the source counts and 
redshift distribution (at the flux limit of the K20 
survey) in the AT-band. This merit function may 
be intended as a modified reduced \ 2 , estimated 
as the sum of the % 2 obtained for each of the four 
observables, and where the number of degrees of 
freedom is replaced by the number of data bins 
(flux or redshift bins, depending on the observ¬ 
able, see Vignali et al.||2009 ). 

The result of our analysis is shown in the top 
panel of Fig.[2j where the merit function is plotted 
as a function of the quantity Zo .5 hi various sep¬ 
arate cases, in which the counts and the redshift 
distributions have been regarded as single observ- 


zini 2006 

). Moreover, 

Ilbert et al. (2013 

) (see also ables, and in the case where these quantities have 

Cimatti et al. 

2006 

Pozzetti et al. 

2010 ), using been considered altogether. 


the new data from the UltraVISTA survey (Mc- 
Cracken et al.|2012 ) in the COSMOS field, showed 


that quiescient galaxies more massive than 10 112 
M 0 do not show any evidence of density evolution 
between 0.8 < 2 < 1.1 and 0.2 < z < 0.5. 

The assumption of z e = 1 is reasonable also in 
the light of the steep decrease at z < 1 of the co¬ 
moving number density of PEP sources associable 


to proto-spheroids, as discussed in Sect. 2.1 


The assumption that star formation in proto¬ 
spheroids is complete at z ~ 1 prevents our model 
to account for the presence of Starbursts, AGN1, 
AGN2 and SF-AGN(SB) objects at z < 1. How¬ 
ever, as visible in Fig. [l] such objects contribute 
by less than ~ 10% to the number density of IR 
sources at 0 < z ^ 1 . 

The quantity z s is the redshift at which star 
formation starts, and is set to z s = 5. We checked 
that z s values larger than 5 do not change our 
results appreciably. 

Our aim is to test whether the proto-spheroid 
formation rate jff- can be described by a piece- 
wise constant function. For the sake of simplicity, 
we considered the most basic piecewise constant 
function, i.e. we divided the redshift range Az 
in two intervals, divided by a value z 0 . 5 , defined 
as the redshift at which half of the proto-spheroid 
population have formed, which is the only free pa¬ 
rameter of our analysis. 

The quantity Z 0.5 is determined by means of a 
‘merit function’, which takes into account the four 


When the redshift distributions and the source 
counts are considered as single observable, the 
minimum of the reduced x 2 falls at Z 0.5 = 3.0 and 
at Z 0.5 = 1.5, respectively. On the other hand, if 
the four observables are considered altogether, our 
analysis indicates a minimum of the merit function 
extending over a region spanning from Z0.5 = 1.5 
and Z0.5 = 3.0. 

In the remainder of this work, we choose for the 
redshift at which half of the proto-spherids have 
form zq _5 ~ 2 as fiducial value, i.e. a value approx¬ 
imately in-between the extremes of the enclosed 
between Z 0.5 = 1.5 and zo .5 = 3.0. 

The assumption of a value for Z0.5 beyond the 
extremes of the distribution causes either a poor 
description of the 160 /jm redshift distribution 
(middle panel of Fig. [2j considering zq .5 = 1.2), or 
a significant underestimation of the 160 /im counts 
at bright fluxes (Sieo^m > 30 mJy, lower panel of 
Fig. [ 2 ] considering z 0 . 5 = 3.5). 

To summarize, we assume that half of the 
proto-spheroids form, at a constant rate, over the 
redshift range 1 < z < Z0.5, and half over the red¬ 
shift interval Z 0.5 < 2 < 5. 

All the results shown in the remainder of this 
work have been obtained with our fiducial assump¬ 
tion, i.e. with Z 0.5 ~ 2. 

3.3. Model results versus observations 

In Fig. [3j the observed differential far-IR 
counts calculated at 160 /rm (upper panel) and the 
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corresponding redshift distribution for the PEP- 
COSMOS data (lower panel) are shown, together 
with the results obtained with our model. 


As shown by the 160 fim band counts, the 
agreement between data and model results is re¬ 
markably good. In particular, the proto-spheroid 
population contributes significantly (being al¬ 
most half of the source counts), in the flux range 
where the differential Euclidean normalized source 
counts peak (Si6o~20 — 30 mJy). 


Also the observed redshift distribution com¬ 
puted at the COSMOS flux limit of 20 mJy (Fig. 
[3} bottom panel) is reproduced by our model. As 
discussed by GPR13, at this flux limit, the com¬ 
pleteness is close to 100 %. 

In Fig. [4j the observed and predicted Ri¬ 
band source counts (upper panel) and redshift dis¬ 
tribution (lower panel) are shown, respectively. 
The model A'-band redshift distribution shown in 
Fig.[4|is compared with the results of the K20 sur¬ 


vey]^ ( Cimatti et al.|[2002| |Mignoli et al.||2005 |. 


The A'-band source counts are well reproduced 
down to A'~23 (Fig. |3J top panel), whereas at 
fainter magnitudes they are underpredicted by our 
model. This result is satisfactory considering that, 
for late-type galaxies, the A'-band counts are cal¬ 
culated from the best-fitting SEDs and from the 
parametric IR LF evolution obtained in GPR13, 
whereas the SED evolution of proto-spheroids is 
computed solely on the basis of the model de¬ 
scribed in Sect. 12.21 


The agreement between our model and data is 
evident also from the computed Af-band source 
redshift distribution (Fig. [4j bottom panel). The 
K20 distribution has been reported since it repre¬ 
sents the most complete spectroscopic survey of a 
AT-band selected sample (i.e. spectroscopic com¬ 
pleteness of ^96 %, Mignoli et al.||2005 ). 

At z>1.5, our analysis indicates that nearly all 
the sources in the redshift distribution are proto¬ 
spheroids. It may be interesting to compare this 
result with the morphological classification of the 
K20 survey (Cassata et al. 2005). Of ~ 300 objects 


2 The K20 survey (area ~52 arcmin 2 ) is composed by 
sources selected in two independent sky regions, one cen¬ 
tred in the Chandra Deep field South (CDFS, ~32.2 
arcmin 2 ) and the second centered around the quasi-stellar 
object 0055-269 (~19.8 arcmin 2 ). More detailed informa¬ 
tion can be found in Cimatti et al. (2002). 


belonging to the CDFS region, ~20-25 sources lie 
at z>1.5, and a large fraction of them (>80 %) 
are classified as irregulars and ellipticals galaxies, 
whereas a small fraction (~ 10 — 20%) are clas¬ 
sified as spiral galaxies. The irregular and ellip¬ 
tical systems may be identified with our proto¬ 
spheroids, caught during their star-forming phase 
and already passive, respectively. Despite the poor 
statistics of the K20 sample at z>1.5, we can re¬ 
gard this result as encouranging. 

In the upper panel of Fig. [5j we show the pre¬ 
dicted differential AT-band counts computed for 
proto-spheroids in various redshift bins. As visible 
also from the redshift distribution plot, the bulk of 
the proto-spheroids is between and z~2, with 
the counts at AT > 19 dominated by spheroids at 
z > 1 and the counts at brighter magnitudes dom¬ 
inated by z <1 sources. Proto-spheroids at higli-^ 

(z > 2) give a significant contribution only at very 
faint magnitudes (A'^23 — 25). As shown in the 
lower panel of Fig. [5j at 160 /./m, instead, most of 
the proto-spheroid population lies at z> 1, peaking 
between 1 < 2 < 2. 

The predicted A'-band and far-IR proto-spheroids 
redshift distribution are significantly different (i.e. 
the A'-band one extending at z< 1 at odd with the 
far-IR), since in our model at 1.6, the star for¬ 
mation is completed in half of the proto-spheroids 
(corresponding to ~0.7 Gyrs after z = 2, see 
Sect. |3.2| ); i.e. at lower redshifts, they will appear 
as red, passive objects. The redshift at which the 
star formation in all spheroids is complete is z^0.8 
(considering the spheroids formed, at the latest 
times, at 2~1). At lower redshifts, the totality of 
these galaxies will evolve as passive systems. 

Later, in Sect. 3.4, we will discuss some im¬ 
plications of our results regarding the stellar mass 
buildup in early-type galaxies. 

3.4. The evolution of the K-band luminos¬ 
ity function 

Fig. [| shows the comparison between the AT- 
band luminosity functions observed at various red¬ 
shifts and the one derived by our model and as¬ 
suming a continuous proto-spheroid formation, as 
discussed in Sect,. 13.21 

In the local Universe, we show the observed to¬ 
tal A'-band LF as determined by |Kochanek et ah] 
|2001[ together with the LF derived for early and 
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late-type galaxies, represented by the red and 
green dashed areas, respectively. At higher red- 
shift, the data are taken from the K20 survey 
up to z ~ 1.75 ( |Pozzetti et al. 2003[ ), and from 
the UKIDSS Ultra Deep Survey (AT~23 over 0.7 
deg 2 , Cirasuolo et al. 2010) in the redshift range 


0.25 < 2 < 3. 


We find an overall good agreement between the 
expected model LF and the data in most of the 
redshift bins of Fig. [G] In the local Universe, 
the early type LF model estimate is consistent 
with the data by construction, since the spheroid 
LF has been normalised to the local A'-band LF. 
On the other hand, the good agreement between 
the computed late-type LF and the data implies 
that the far-IR selection is not missing a signifi¬ 
cant number of sources classified as late-types also 
in the NIR, and that the SED-fitting procedure 
used for the classification is nearly accurate. The 
slightly shallower decline, at bright A'-band mag¬ 
nitudes, of the computed late-type LF with re¬ 
spect to the observed one reflects the shape of IR- 
derived luminosity functions, typically well fitted 
by a modified Schechter function rather than by 


a pure Schechter (e.g. Saunders et al. 1990 Pozzi 

etaU|2004 

Salim & Lee 2012 

)• 


We find that the late-type population exhibits a 
steeper faint-end slope than early-types, and that 
they dominate the LF at the fainter luminosities 
up to 2 ~ 2. 

The largest discrepancy between the model pre¬ 
dictions and the data is visible in the redshift 
range 1.75 tzt 2.75, where our model overpredicts 
the number of luminous galaxies, at the bright end 
of the LF {M k < - 24.5). 

Possible reasons for the discrepancy are the sim¬ 
plicity of our model (i.e. we recall that we are us¬ 
ing only one proto-spheroid model for the whole 
luminosity range, see Sect. 2.2) and the cosmic 
variance, whose effects are visible also from the dif¬ 
ference between the bright-end LF of the K20 sur- 


vey ( 

Pozzetti et al.|2003 

) and that of the UKIDSS 

survey ( 

Cirasuolo et al. 

2010 

)■ 


At high redshift (2.7 <; z <; 3.25), the agree¬ 
ment between model results and the observed LF 
is remarkable. 


Considering the whole redshift range, our re¬ 
sults indicate that the contribution of late type 
galaxies to the total LF is significant up to 2~2. 





Fig. 7. - Comparison between the observed and 
predicted rest-frame SEDs for proto-spheroids in 
the COSMOS field, calculated at the 160 (im sur¬ 
vey flux limit, and at three different redshifts 
(from top to bottom: z = 1.5, z = 2.0 and z = 
2.5). In each panel, the observed SEDs are shown 
as crosses and represents the sources selected in 
the redshift range 1.3 < z < 1.7, 1.8 < 2 < 2.2 
and 2.3 < 2 < 2.7. The model SEDs have been 
calculated as described in Sect. 3.5 and are repre¬ 
sented by the colour-coded regions, which express 
the expected number of sources at wavelength A 
and with flux S , normalized to the total number 
of sources at that wavelength . The green areas 
represent regions where the probability of finding 
a source is greater than 10 %. 
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At 2 > 2, the A'-band luminosity function is dom¬ 
inated by proto-spheroids. 

3.5. Observed and predicted SEDs 



Fig. 8 .— The observed local early-type A-band 
LF from Kochanek et al. (2001). The red short- 
dashed line and the red shaded area represent the 
Schechter fit to the local A"-band LF and the as¬ 
sociated l -<7 uncertainties respectively. The long- 
dashed black vertical lines mark the portions of 
the local LF associated to the three mass models, 
as described in Sect. 13.61 


As further check, in this Section we aim at test¬ 
ing how the SEDs of our proto-spheroid model 
galaxies compare with the observed SEDs of the 
non-spiral galaxies of the GPR13 sample. 


In Fig. [7J the theoretical proto-spheroid rest- 
frame SEDs are compared with the observed SEDs 
at high redshift (z>1.5). The comparison has 
been performed in the COSMOS field, consider¬ 
ing the 160 (im flux limit of 20 mJy, at which 
the completeness of the survey is close to 100 % 
(see Sect. 3.3). The model predictions have been 
computed at three redshifts (z=1.5, 2 and 2.5), 
while the observed SEDs have been selected in 
the redshift bins 1.3< 2 <1.7, 1.8< 2 <2.2 and 
2.3< 2 <2.7. 


The expected model SEDs have been simu¬ 
lated, at each redshift, starting from the theo¬ 
retical proto-spheroids luminosity function calcu¬ 
lated at 160 pm. By construction, the total proto¬ 
spheroids LF at a given redshift 2 is the sum of 
LF of the proto-spheroid sub-populations formed 
at different epochs and present at that redshift. At 
each luminosity and at each redshift, we can define 
a weight Wi relative to an extant sub-population, 
given by the ratio between its LF and the total 
LF at that redshift. The SEDs plotted in Fig. [7] 
are those of the extant i —th sub-populations at 
2 = 1.5, 2 = 2 and 2 = 2.5, each multiplied by its 
weight Wi. 

Overall, the theoretical SEDs are in satisfac¬ 
tory agreement with the observations of GPR13, 
in particular at 2 = 2.5. At lower redshift, the 
model SEDs are a factor ~ 5 larger than the ob¬ 
served ones in the Mid-IR regime, indicating a 
warm-dust temperature higher that that inferred 
from the observed galaxies. However, this discrep¬ 
ancy does not affect the main results of this study, 
which is based on Near-IR and Far-IR data, where 
the observed SEDs are in good agreement with the 
model results. 


3.6. The effects of downsizing 
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The single-mass model described above is de¬ 
signed to describe the evolution of the bulk of 
the population of spheroids. A more sophisticated 
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Fig. 9.- Observed and predicted iO-band rest-frame luminosity function in eight different redshift bins. 
The solid lines represent the total A"-band LFs calculated taking into accounts the effects of downsizing of 
proto-spheroids, as described in Sect. |3.6] The dashed lines represent the contribution of late-type galazies. 
The shaded areas, solid circles and open squares are as described in Fig. [G] 
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Fig. 10.— Upper panel: differential extragalac- 
tic source counts normalized to the Euclidean 
slope at 160 /im. The curves represent the model 
predictions, calculated taking into account the 
downsizing of proto-spheroids. The thin solid 
line with crosses, the thick solid line and the 
dot-dashed line represent the contribution of the 
proto-spheroid population of present-day stellar 
mass 3 x 10 10 Mq, 10 10 Mq and 10 12 Mq, respec¬ 
tively; the dotted green line represents the spi¬ 
ral, the dot-dot-dot-dashed dark-green the SF- 
AGN(Spiral) while the long-dashed dark-green the 
spiral+SF-AGN(Spiral) populations. The obser¬ 
vational data are as described in the top panel of 
Fig. [3] In the lower panel, the redshift distribu¬ 
tion is shown. The curves are as in the top panel 
and the observational data are as described in the 
bottom panel of Fig. [3] 



z 


Fig. 11.— Upper panel: AT-band differential ex- 
tragalactic source counts calculated taking into ac¬ 
count the effects of downsizing in proto-spheroids 
as described in Sect. 13.61 The curves are as in 
Fig. [TO} The observational data as described in 
Fig. Lower panel: Predicted redshift distri¬ 
bution of the AT-band selected sources calculated 
taking into account the effects of downsizing. The 
curves are as above and the observational data are 
as described in Fig. [I] 
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treatment should include proto-spheroid models 
of various masses, suited to study the differen¬ 
tial evolution of populations of various masses, in¬ 
cluding their downsizing behaviour. To this pur¬ 
pose, we use chemo-spectrophotometric models of 
proto-spheroids of three different present-day stel¬ 
lar masses, and we investigate the effects of down¬ 
sizing on the main observables considered in this 
paper. 


The three models are characterised by present- 
day stellar masses of 3 x 10 1 o Mq, 10 11 Mq and 
10 12 Mq. The main parameters describing the 
models are the star formation efficiency (from 5 
up to 20 Gyrs -1 for masses between 3 x 10 10 and 
10 12 Mq) and the time of the onset of the galactic 
wind (from 0.3 Gyrs up to 0.9 Gyrs). For fur¬ 
ther details, the reader is referred to [Calura et al. 


(2014), where the models are fully described. 

In analogy with the method described in 
Sect. |3.2[ the present-day K —band luminosity 
function of Kochanek et al. (20011 is used to nor¬ 
malize the total number density of local spheroids. 
In this case, the present-day early-type AT-band 
LF is divided into three parts, each one corre¬ 
sponding to a different population. Each pop¬ 
ulation is then evolved backwards by means of 
the corresponding SED, as earlier computed with 
GRASIL, starting from the star formation his¬ 
tory and from the chemical evolution. The three 
models of different masses are associated to the 
three proto-spheroid populations as follows. The 
evolution of the faintest galaxies, i.e. those at 
AT—band magnitudes M k[AB) ^ 21, are described 
by the 3 x 10 10 Mq model. The intermediate 
mass model, characterised by a present-day stel¬ 
lar mass of 10 11 M 0 , describes the evolution of 
galaxies with magnitudes 21 xi M k[AB] ^ 23. 
Finally, the evolution of the brightest galaxies 
(23 Xi M k[AB]) is described by the highest-mass 
model, with present-day stellar mass of 1O 12 M 0 . 
The stellar masses have been converted to K— 
band magnitudes by means of the mass-to-light 
ratios (t). calculated for each model as a func¬ 
tion of the age. 


In this case, for each population, the functional 
form of the formation rate is not assumed a pri¬ 
ori; this approach would generate a considerably 
wide parameter space and a significant degeneracy 
in the results. 


In the last few years, various observational in¬ 


vestigations of early-type galaxies have been fo¬ 
cused on the determinations of the ages of their 
stellar populations (Renzini]20061, providing clear 
relations between present-day stellar masses and 
formation redshift. Here we assume the relation 
between formation epoch and mass as determined 
by Moresco et al. (2010), based on the multi¬ 
band SED-fitting analysis of a zCOSMOS sample 
of early type galaxies up to z = 1. The work of 


Moresco et al. (2010) present measures of the look¬ 


back time to formation for galaxies of stellar mass 
between 1O 1O M 0 and 1O 11 M 0 . 

In our picture, for each population the forma¬ 
tion rate is given by a gaussian function cen¬ 
tered at the age measured by Moresco et al. ( 201 0), 
where the standard deviation corresponds to the 
associated age dispersion. It is worth stressing 
that the ages measured by Moresco et al. (2010) 
have been converted to our cosmology. 

For galaxies of 1O 12 M 0 , we assume a mean 
present-day age of 12 Gyrs, in agreement with 
other independent studies (see, e. g., Renzini 
2006), and we adopt the same 1 — a uncertainty 
as the one of the systems of mass IO^Mq. 


The three AT-band magnitudes used to divide 
the local A'-band LF have been calculated from 
the present-day SEDs of our models, computed at 
the mean age of each population, i.e., following 
Moresco et al. (20101, at ~ 7.8 Gyr, 10.5 Gyr and 


12 Gyrs for the 3 x 1O 1O M 0 , 1O 11 M 0 and 1O 12 M 0 
model, respectively. 


In Fig. [8j we show our division of the local 
AT-band LF. The vertical dashed lines mark the 
extremes of the three intervals the LF has been 
divided into. Each line corresponds to a magni¬ 
tude value falling half-way between the present- 
day magnitudes of two contiguous models. The 
total population of spheroids has then been di¬ 
vided into three populations of different masses. 
In analogy with our single-mass model, the back¬ 
wards evolution of each population is determined 
by the evolution of the SED of the corresponding 
model. 


The evolution of the A'-band LF obtained by 
taking into account the effects of downsizing is 
shown in Fig. [9] At redshift z > 2, the observed 
LF is now slightly underestimated at the break 
magnitudes, i.e. at Mk ~ —22. In principle, this 
discrepancy could be reduced by fine-tuning the 
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adopted mean age of the 3 x 10 iu and 10 ii M Q 
proto-spheroids, e. g. by slightly increasing this 
quantity, in agreement with the results of other 


studies of local early type galaxies (see Renzini 


2006). However, we choose to use the age-mass re¬ 


lation of Moresco et al. (2010) at face value and to 
not regard this disagreement as particularly wor¬ 
rying, considering the uncertainty of the present- 
day age of such objects, of the order of f Gyr 
( Moresco et al.|[2010 ). 


Another remarkable effect of downsizing in the 
predicted high -2 LF is in the slope of the faint 
end: in our picture, the least massive progenitors 
are very rare at z > 2 and this causes a strong 
fall-off of the LF at Mk > —22. 


In general, at lower redshift, the predicted total 
A'—band LF is in good agreement with the obser¬ 
vations. 


The effects of the proto-spheroids of various 
masses are visible also in Fig. [TOJ where we show 
the far-IR differential counts and redshift distribu¬ 
tion. From the top panel Fig. [lOj we can see that 
the total proto-spheroidal contribution is similar 
to the ones computed with the single-mass model, 
shown in Fig. [3] 


The 10 n M Q mass model gives a dominant con¬ 
tribution to the far-IR source counts, in particular 
in the flux range between ~ 2 mJy and ~ 30 mJy, 
and to the redshift distribution from z ~ 0.8 to 
2 ~ 2 (lower panel of Fig. 10). 

It is also worth noting that the introduction of 
downsizing allows us to reproduce slightly better 
the observed redshift distribution of far-IR sources 
at z > 2.5, thanks to the contribution of the most 
massive systems (i.e. 1O 12 M 0 ). 


The predicted AT—band differential counts and 
the redshift distribution of K-band sources are 
shown in the top and bottom panel of Fig. El 
respectively. In general, the total counts are very 
similar to the ones shown in Fig. [4] (upper panel, 
i. e. assuming a single-mass model to describe the 
whole population of proto-spheroids), with the ex¬ 
ception of the faint-end, where the disagreement 
between data and model results is slightly more 
marked. Once again, this effect, could be reduced 
by fine-tuning the mean age of the 3 x 10 10 and 
lO n M 0 mass model (i.e. moving the peak of for¬ 
mation at higher redshift). 

The similarity of the results computed with the 


single-mass model and obtained with the three- 
mass models are in agreement with other previ¬ 
ous chemical evolution studies of early-type galax¬ 
ies. In fact, Matteucci & Vettolani (1988) showed 
that the chemical enrichment of the intra-cluster 
medium is dominated by galaxies with L ~ A*. 

The 10 11 M 0 mass model dominates the K- 
band proto-splreroids source counts from the 
brighter magnitudes down to A'^22, while at 
fainter magnitudes the low mass populations (i.e. 
3 x 1O 1O M 0 ) becomes the prevailing one. 

Consistently with the source counts, the K- 
band redshift distribution (at the limit magni¬ 
tude of A'~21.9) shows that the proto-spheroids 
of 10 11 M 0 give the main contribution at z > 0.5; 
at lower redshifts, the major contribution is given 
by the lowest mass population (i.e. 3 x 1O 1O M 0 ). 
The highest mass model (i.e. 1O 12 M 0 ) population 
is always negligible (given the small area of the 
K20 survey and their low volume density); in the 
source counts, this population gives a contribution 
only at very bright magnitudes (K < 19). 

It is worth stressing that the formation rate pro¬ 
vided by the age measures of Moresco et al. (2010) 
allow us to reproduce the observed distributions 
with good accuracy: no fine-tuning of any forma¬ 
tion rate was required in this case. 


Our main conclusion is that at most magni¬ 
tudes/fluxes, the main contributors of the total 
differential far-IR and A'—band counts of proto¬ 
spheroids are represented by a population of 
present-day mass which corresponds roughly to 
the break of the local early-type stellar mass func¬ 
tion. 


3.7. Formation redshift distribution and 
mass assembly 

In Sections 13.21 and 13.61 we have shown that it 
is possible to reproduce a set of multi-wavelength 
galactic observables across a wide redshift range 
(0 < z < 3) by means of a population of late-type 
galaxies, whose behaviour is determined empiri¬ 
cally by the evolution of their far-IR-LF, plus a 
population of proto-spheroids. In Sect. |3.2[ we 
have assumed a piecewise constant formation rate 
across the redshift range 1 < 2 < 5, whereas in 
Sect. |3.6| we have divided the proto-spheroids in 
three classes of different mass and, for each class, 
we have assumed an observationally-derived age- 
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mass relation. In Fig. [l2j we show the evolution 
of the formation rate dN/dzf orm computed in the 
case of the single-mass model (dashed line), to¬ 
gether with the formation rates of the three popu¬ 
lations of spheroids described in Sect. 3.6 (dashed 
areas). These latter curves have been normal¬ 
ized according to the local stellar mass density 
of the corresponding populations and the sum of 
their integrals has been normalized to one as for 
the single-mass model. We find that the com¬ 
bined contributions of low-mass (3xl0 10 M 0 ) and 
high-mass (10 12 M 0 ) spheroid populations sum up 
to ~35% of the present stellar mass density (we 
adopted 10 8 to 10 13 M 0 as extremes for the inte¬ 
gration). 

As visible in the plot, the dN/dzf orm function 
calculated in the case of the single-mass model is 
in fairly good agreement with the one correspond¬ 
ing to M* = 10 11 Mq of the three-mass model (i.e. 
in both cases, at 2 half of the total present- 
day mass has already formed) and obtained by a 
completely independent method. The global con¬ 
sistency of the two curves, together with the dom¬ 
inance of the M* = 1O 11 M 0 population at z < 2.5 
(see Sect. 3.6), outlines the validity of our basic 
assumption, i.e. that the evolution of the whole 
proto-spheroid population can be traced by the 
evolution of galaxies at the break of the present- 
day early type AT-band LF. 

The most direct observables suited to study the 
mass assembly history of spheroids is the evolution 
of their stellar mass density (SMD, i.e. Dickinson 


et al. 2003). From the theoretical point of view, 


studying how the stellar mass of galaxies grows 
with time is fundamental to assess the role of var¬ 
ious processes in determining one of the most basic 
galactic parameters, including the roles of mergers 
and the relative roles of various quenching mech¬ 
anisms in driving the galactic star formation his¬ 
tories. 

On the other hand, one should bear in mind 
that the stellar mass in not a direct observable 
as its observational determination depends on a 
number of model-dependent assumptions, includ¬ 
ing the a-priori choice of an initial mass function, 
a star formation history and a metallicity: these 
factors constitute major sources of uncertainty in 
the observational determination of the stellar mass 
in galaxies (see e.g. Pozzetti et al.|2007 Maraston 
etaLpOLOl ). 


To calculate this quantity, we convert the 
A'-band LFs as calculated at various redshifts 
(Sect. [3)4| and shown in Fig. [6] into stellar mass 
functions (MFs) by means of the mass-to-light ra¬ 
tios it) computed at the time t and, at each 
epoch, by properly taking into account the contri¬ 
butions of galaxies born at different epochs. 


Our mass-to-light ratios range from R ~ 
0 .2Mq/Lq^k, characterising the end of the SF 
epoch, up to a present-day value of R ~ 1Mq/Lq } k 
Such values are comparable to those used by other 
authors (see e.g. Drory et al.||2004 Narayanan & 


Dave 2013) in theoretical studies of the stellar 


populations of early type galaxies. 


The SMD is calculated from the integral of the 
stellar mass function 0(M+, z) over the mass range 
10 s to 10 13 M 0 , according to: 


r 13 

P*(z)= / </>(M', z)M' if d logMl . (6) 

J 8 


In the bottom panel of Fig. |13[ we compare 
the evolution of the stellar mass density calculated 
for spheroids with various estimates from the lit¬ 
erature (see caption of Fig. 13 for details). The 


SMD from Ilbert et al. (2013) has been derived in 


the COSMOS field from a AT-band selected sam¬ 
ple down to K=2A and using the new UltraVISTA 
DR1 data release. In Fig. 0 we show the total 
SMD as derived by Ilbert et al. (2013), i.e. in¬ 
cluding both quiescent and late-type systems; the 
original values have been divided by 0.55 to con¬ 


vert them from a Chabrier (2003) to a Salpeter 
(1955) initial mass function. In Fig. we also 


show a compilation of observational data by |San- 


tini et al. (2012). 


In our scenario, at z > 2 the cosmic star for¬ 
mation is dominated by proto-spheroids; in this 
redshift range, the SMD calculated in this work is 
consistent with the total estimates of lllbert, et al.l 
(2013) and Santini et al. (20121 within the uncer¬ 
tainties. This result supports a scenario in which 
most of the mass at high-redshift (i.e. observed at 
z > 2) is in proto-spheroids. Moreover our results 
indicate that half of the proto-spheroids mass must 
have formed at z > 2, and the remaining between 
1 and z~2. 

At lower redshift, the SMD values of |Ilbert et al.| 
(2013) and Santini et al. (20121 are slightly larger 
than our estimates; this result is not surprising, 
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as our estimates are lacking the contribution from 
late-type galaxies. 

A detailed study of the contribution of late-type 
systems to the SMD evolution requires the con¬ 
struction of a library of template SEDs similar to 
those used in this paper to study proto-spheroids; 
such a task is deferred to future work. 


4. Conclusion 


In this paper we have interpreted a set of multi¬ 
wavelength observations of galaxies across a large 
fraction of the cosmic epoch (0 < z < 3) by means 
of a new phenomenological approach. Our ap¬ 
proach matches a ‘backward’ parametric evolution 
of late-type galaxies based on the observed IR lu¬ 
minosity functions of GPR13 to a model for proto¬ 
spheroids where the spectro-photometric evolution 
of dust has been calculated self-consistently, on 
the basis of idealised star formation histories able 
to account for a large set of observational data. 
These data include the dust budget observed in lo¬ 
cal passive galaxies and detected in high-redshift 


massive starbursts (CPM08; Pipino et al. 2011) 


and their SEDs as observed in both the local and 
distant universe ( |Schurer et al. 2009). This ap¬ 
proach has been developed in order to simulta¬ 
neously account for observables related mostly to 
the star-formation process and young stellar pop¬ 
ulations, traced by far-IR emission, and also for 
the ones associated to evolved stellar populations, 
whose presence is detected in the near-IR. 

As a first step, following GPR13, IR sources have 
been divided into late-type and potential progeni¬ 
tors of present-day spheroids on the basis of their 
SED and of the evolution of their comoving num¬ 
ber density. The number density of late-type 
galaxies has been parameterized as in GPR13. 
The number density of proto-spheroids has then 
been normalized according to the local A'-band 
LF of early-type galaxies from Kochanek et al.| 
(2001). Finally, the evolution of their A'-band 
and far-IR LFs has been calculated by means 
of the model SEDs from the physical chemo- 
spectrophotometric model described in Schurcr 
et al. (20091. We have also tested the effects of 


mass downsining on the observables studied in 
this paper, by dividing the proto-spheroid popu¬ 
lation into three populations of different masses, 


and tracing their backwards evolution by means 
of chemo-spectrophotometric models of galaxies of 
present-day stellar mass of 3 x 1O 1O M 0 , lO n M 0 
and lO 12 Af 0 . Our main conclusions can be sum¬ 
marised as follows. 

• We have considered a continuous formation 
for proto-spheroids, occurring across the red- 
shift range 1 < z < 5. In this picture, the 
proto-spheroid formation is described by a 
piecewise-constant function, where the key 
parameter is the redshift Z 0.5 at which half 
the population has formed. To constrain this 
parameter, we have performed a statistical 
test by taking into account the following four 
observables: the source counts and redshift 
distribution (at the flux limit of the COS¬ 
MOS field) at 160 pm and the source counts 
and redshift distribution (at the flux limit of 
the K20 survey) in the AT-band. Our anal¬ 
ysis shows that our full set of observables 
can be reproduced by assuming zo .5 in the 
interval 1.5 < zo .5 < 3. The adoption of val¬ 
ues for zq .5 lower than 1.5 implies a poor 
description of the 160 pm redshift distribu¬ 
tion; on the other hand, the assumption of 
zq .5 larger than 3 causes a poor description 
of the 160 pm counts at the brightest fluxes 
(S 16 o^m > 30 mJy). 

By assuming zq.$ ~ 2 as fiducial value, 
we find that at the far-IR flux limit of the 
PEP-COSMOS survey, all the PEP-selected 
sources at z > 2 can be explained as progen¬ 
itors of local spheroids caught during their 
formation. 

• A continuous formation of proto-spheroids 
with zo .5 ~ 2 allows us to reproduce the 
differential extragalactic source counts ob¬ 
served at 160 p?n, with proto-spheroids 
contributing significantly at the flux range 
where the differential Euclidean normalized 
source counts peak (<5i6o~20—30 mJy). Also 
the differential counts in A'-band are repro¬ 
duced by our scenario, with the exception 
of the faint-end portion. This discrepancy 
is probably due to the non-evolving slope 
of the proto-spheroid luminosity function as 
adopted in this work. 
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• The differential AT-band counts computed in 
various redshift bins indicate that the bulk 
of the proto-spheroids he between z ~0 and 
z~2, with the counts at K > 19 dominated 
by spheroids at z > 1 and the counts at 
brighter magnitudes dominated by z < 1 
sources. High-z proto-spheroid sources con¬ 
tribute significantly to the AT-band counts 
only at very faint magnitudes (AW23 — 25). 
From the far-IR counts, we have seen that 
most of the star-forming proto-spheroid pop¬ 
ulation lies between 1 < z < 2 . As further 
check, we found a satisfactory agreement 
between the simulated SEDs of the model 
proto-spheroids and the spectra of non-spiral 
sources in the COSMOS field at z > 1.5. 

• Our results obtained by taking into account 
the effects of downsizing show that the main 
contribution to the counts comes from ob¬ 
jects of present-day stellar mass of 10 11 A/ Q , 
corresponding roughly to the break of the 
local early-type stellar mass function. 

• Finally, we have studied the evolution of 
the stellar mass density as implied by our 
results and we have compared it to other 
estimates from the literature. At high-z 
(z > 2 ), we find a good consistency between 
the SMD calculated in this work and the to¬ 
tal estimates from other authors, support¬ 
ing a scenario in which most of the mass at 
high-redshift is in proto-spheroids. More¬ 
over our results indicate that half of the 
proto-spheroids mass must have formed at 
z > 2, and the remaining ones between z~l 
and z~2. The same results are obtained also 
when downsizing is taken into account. This 
confirms the validity of our approach, i.e. 
that the mass growth of the proto-spheroid 
population can be modelled by means of a 
single-mass model, describing the evolution 
of galaxies at the break of the present-day 
early-type A'-band LF. 

At lower redshift, the total SMD estimates 
from literature are slightly larger than our 
estimates. This is likely due to the fact that 
late-type galaxies have not been considered 
in our study. In the future, the contribution 
of late-type galaxies to the SMD will need 
to be assessed in order to provide more accu¬ 


rate predictions on the evolution of the SMD 
and to gain further hints on their formation 
history. 


As a final note, we stress that our results do not 
allow us to quantitatively assess the role of merg¬ 
ers as main drivers of the star formation history 
of proto-spheroids, which is probably dominant 


at large redshifts (Renzini 2009 Conselice et al. 


2009). A merger-driven formation of spheroids is 


still in agreement with our results, but if the as¬ 
sembly of the sytems of stellar mass 1 O 11 M 0 oc¬ 
curs preferentially at z >1. Presently, the major 
difficulty of cosmological scenarios, where proto¬ 
spheroids form as the results of a sequence of merg¬ 
ing episodes, is to have the starbursts associated 
to the most massive systems completing their star 
formation histories before ~ 1 Gyr, a duration 
which allows to account for a wide set of multi¬ 
wavelength properties at high redshift as seen here 
and in other works (see, e.g., Granato et al.|[2004 


Silva et al.|2 005, Lapi et al. 2011 ) , as well as their 

main scaling relations as traced by the basic stel¬ 


lar population diagnostics of local ellipticals (Mat- 
teucci|[T994 Renzini|[2006 ). 

To study the IR properties of galaxies by 
means of cosmological models, one needs to take 
into account the dust production processes self- 
consistently as in the model used in this work. 
To accomplish such a task, the next-generation 
of A—CDM semi-analytical models will have to 
include a detailed treatment of the chemical evo¬ 
lution of refractory elements (Calura & Mend 


2011 ). 
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Fig. 12.— Differential proto-spheroid formation 
rate computed in the case of the single-mass model 
(dashed line), together with the formation rates of 
the three populations (dashed areas). Each curve 
is normalized to the present-day stellar mass den¬ 
sity of the corresponding population. The curve of 
the single mass model is normalized to the total 
present-day stellar mass in early-type galaxies. 


z 



Age of the Universe (Gyr) 


Fig. 13.— Observed and predicted evolution of 
stellar mass density as a function of cosmic time. 
The solid and dashed lines are the theoretical 
SMD calculated for the one-mass and the three- 
mass models, respectively. The grey circles are 
the SMD values from the compilation of |Santim| 


et al. (2012|. The dashed region represents the 


estimates from Santini et al. (2012) that take into 
account all the glocal sytematic uncertainties. The 
filled grey triangles represent the total SMD as de¬ 
rived in Ilbert et al. (2013). 


23 










































REFERENCES 

Baldry, I. K., Glazebrook, K., Brinkmann, J., 
et al. 2004, ApJ, 600, 681 

Berta, S., Magnelli, B., Nordon, R., et al. 2011, 
A&A, 532, A49 

Bethermin, M., Dole, H., Beelen, A., & Aussel, H. 
2010, A&A, 512, A78 

Betlreraiin, M., Dole, H., Lagache, G., Le Borgne, 
D., & Penin, A. 2011, A&A, 529, A4 

Bonato, M., Negrello, M., Cai, Z.-Y., et al. 2014, 
MNRAS, 438, 2547 

Burgarella, D., Buat, V., Gruppioni, C., et al. 
2013, A&A, 554, A70 

Cai, Z.-Y., Lapi, A., Xia, J.-Q., et al. 2013, ApJ, 
768, 21 

Calura, F., Gilli, R., Vignali, C., et al. 2014, MN¬ 
RAS, 438, 2765 

Calura, F., & Mend, N. 2011, MNRAS, 413, LI 

Calura, F., Pipino, A., & Matteucci, F. 2008, 
A&A, 479, 669 (CPM08) 

Caputi, K. I., Lagache, G., Yan, L., et al. 2007, 
ApJ, 660, 97 

Cassata, P., Cimatti, A., Franceschini, A., et al. 
2005, MNRAS, 357, 903 

Chabrier, G. 2003, ApJL, 586, L133 

Cimatti, A., Daddi, E., & Renzini, A. 2006, A&A, 
453, L29 

Cimatti, A., Mignoli, M., Daddi, E., et al. 2002, 
A&A, 392, 395 

Ciotti, L. 2009, Nuovo Cimento Rivista Serie, 32, 
1 

Cirasuolo, M., McLure, R. J., Dunlop, J. S., et al. 
2010, MNRAS, 401, 1166 

Conselice, C. J., Yang, C., & Bluck, A. F. L. 2009, 
MNRAS, 394, 1956 

Cook, M., Lapi, A., & Granato, G. L. 2009, MN¬ 
RAS, 397, 534 


Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, 
J. G. 1996, AJ, 112, 839 

Cresci, G., Hicks, E. K. S., Genzel, R., et al. 2009, 
ApJ, 697, 115 

Cucciati, O., Tresse, L., Ilbert, O., et al. 2012, 
A&A, 539, A31 

De Lucia, G., Fontanot, F., Wilnran, D., & 
Monaco, P. 2011, MNRAS, 414, 1439 

Dekel, A., Birnboim, Y., Engel, G., et al. 2009, 
Nat, 457, 451 

Dickinson, M., Papovich, C., Ferguson, H. C., & 
Budavari, T. 2003, ApJ, 587, 25 

Dole, H., Rieke, G. H., Lagache, G., et al. 2004, 
ApJS, 154, 93 

Draine, B. T., & Li, A. 2007, ApJ, 657, 810 

Driver, S. P., Robotham, A. S. G., Bland- 
Hawthorn, J., et al. 2013* MNRAS, 430, 2622 

Drory, N., Bender, R., Feulner, G., et al. 2004, 
ApJ, 608, 742 

Elbaz, D., Cesarsky, C. J., Fadda, D., et al. 1999, 
A&A, 351, L37 

Fontanot, F., & Somerville, R. S. 2011, MNRAS, 
416, 2962 

Fontanot, F., Somerville, R. S., Silva, L., Monaco, 
P., & Skibba, R. 2009, MNRAS, 392, 553 

Forster Sclrreiber, N. M., Genzel, R., Bouche, N., 
et al. 2009, ApJ, 706, 1364 

Franceschini, A., Rodighiero, G., Vaccari, M., 
et al. 2010, A&A, 517, A74 

Frayer, D. T., Fadda, D., Yan, L., et al. 2006, AJ, 
131, 250 

Gargiulo, I. D., Cora, S. A., Padilla, N. D., et al. 
2014, ArXiv e-prints, arXiv:1402.3296 

Genzel, R., Burkert, A., Bouche, N., et al. 2008, 
ApJ, 687, 59 

Gnerucci, A., Marconi, A., Cresci, G., et al. 2011, 
A&A, 528, A88 

Granato, G. L., De Zotti, G., Silva, L., Bressan, 
A., & Danese, L. 2004, ApJ, 600, 580 


24 



Gruppioni, C., Lari, C., Pozzi, F., et al. 2002, MN- 
RAS, 335, 831 

Gruppioni, C., Pozzi, F., Zamorani, G., & Vignali, 
C. 2011, MNRAS, 416, 70 

Gruppioni, C., Pozzi, F., Andreani, P., et al. 2010, 
A&A, 518, L27+ 

Gruppioni, C., Pozzi, F., Rodighiero, G., et al. 
2013, MNRAS, 432, 23 (GPR13) 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 
142 

Ilbert, O., McCracken, H. J., Le Fevre, O., et al. 
2013, A&A, 556, A55 

Kochanek, C. S., Pahre, M. A., Falco, E. E., et al. 
2001, ApJ, 560, 566 

Kodama, T., Yamada, T., Akiyama, M., et al. 
2004, MNRAS, 350, 1005 

Lacey, C. G., Baugh, C. M., Frenk, C. S., et al. 
2008, MNRAS, 385, 1155 

Lapi, A., Gonzalez-Nuevo, J., Fan, L., et al. 2011, 
ApJ, 742, 24 

Lara-Lopez, M. A., Cepa, J., Bongiovanni, A., 
et al. 2010, A&A, 521, L53 

Larson, R, B. 1976, MNRAS, 176, 31 

Lutz, D., Poglitsch, A., Altieri, B., et al. 2011, 
A&A, 532, A90 

Madau, P., & Pozzetti, L. 2000, MNRAS, 312, L9 

Magnelli, B., Elbaz, D., Chary, R. R., et al. 2009, 
A&A, 496, 57 

Magnelli, B., Popesso, P., Berta, S., et al. 2013, 
A&A, 553, A132 

Mannucci, F., Cresci, G., Maiolino, R., Marconi, 
A., & Gnerucci, A. 2010, MNRAS, 408, 2115 

Maraston, C., Pforr, J., Renzini, A., et al. 2010, 
MNRAS, 407, 830 

Marchesini, D., van Dokkum, P. G., Forster 
Schreiber, N. M., et al. 2009, ApJ, 701, 1765 

Matteucci, F. 1994, A&A, 288, 57 


Matteucci, F., & Vettolani, G. 1988, A&A, 202, 
21 

McCracken, H. J., Milvang-Jensen, B., Dunlop, J., 
et al. 2012, A&A, 544, A156 

Mignoli, M., Cimatti, A., Zamorani, G., et al. 
2005, A&A, 437, 883 

Moresco, M., Pozzetti, L., Cimatti, A., et al. 2010, 
A&A, 524, A67 

Narayanan, D., & Dave, R. 2013, MNRAS, 436, 
2892 

Niemi, S.-M., Somerville, R. S., Ferguson, H. C., 
et al. 2012, MNRAS, 421, 1539 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Oliver, S. J., Wang, L., Smith, A. J., et al. 2010, 
A&A, 518, L21+ 

Papovich, C., Dole, H., Egami, E., et al. 2004, 
ApJS, 154, 70 

Pilbratt, G. L., Riedinger, J. II., Passvogel, T., 
et al. 2010, A&A, 518, LI 

Pipino, A., Fan, X. L., Matteucci, F., et al. 2011, 
A&A, 525, A61 

Pipino, A., Kawata, D., Gibson, B. K., & Mat¬ 
teucci, F. 2005, A&A, 434, 553 

Pipino, A., Matteucci, F., Borgani, S., & Biviano, 
A. 2002, New A, 7, 227 

Poggianti, B. M. 1997, A&A Supp., 122, 399 

Polletta, M., Tajer, M., Maraschi, L., et al. 2007, 
ApJ, 663, 81 

Porter, L. A., Somerville, R. S., Primack, J. R., 
& Johansson, P. H. 2014, ArXiv e-prints, 
arXiv: 1407.0594 

Pozzetti, L., Cimatti, A., Zamorani, G., et al. 
2003, A&A, 402, 837 

Pozzetti, L., Bolzonella, M., Lamareille, F., et al. 
2007, A&A, 474, 443 

Pozzetti, L., Bolzonella, M., Zucca, E., et al. 2010, 
A&A, 523, A13 

Pozzi, F., Gruppioni, C., Oliver, S., et al. 2004, 
ApJ, 609, 122 


25 



Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Puget, J.-L., Abergel, A., Bernard, J.-P., et al. 
1996, A&A, 308, L5+ 

Renzini, A. 2006, ARA&A, 44, 141 
2009, MNRAS, 398, L58 

Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., 
et al. 2009, ApJ, 692, 556 

Rodighiero, G., Vaccari, M., Franceschini, A., 
et al. 2010, A&A, 515, A8 

Rodighiero, G., Daddi, E., Baronchelli, I., et al. 

2011, ApJL, 739, L40 

Sajina, A., Yan, L., Fadda, D., Dasyra, K., & 
Huynh, M. 2012, ApJ, 757, 13 

Salim, S., & Lee, J. C. 2012, ApJ, 758, 134 

Salpeter, E. E. 1955, ApJ, 121, 161 

Santini, P., Fontana, A., Grazian, A., et al. 2012, 
A&A, 538, A33 

Saunders, W., Rowan-Robinson, M., Lawrence, 
A., et al. 1990, MNRAS, 242, 318 

Schawinski, K., Simmons, B. D., Urry, C. M., 
Treister, E., & Glikman, E. 2012, MNRAS, 425, 
L61 

Schmidt, M. 1959, ApJ, 129, 243 

Schurer, A., Calura, F., Silva, L., et al. 2009, MN¬ 
RAS, 394, 2001 

Shimasaku, K., & Fukugita, M. 1998, ApJ, 501, 
578 

Shupe, D. L., Rowan-Robinson, M., Lonsdale, 
C. J., et al. 2008, AJ, 135, 1050 

Silva, L., De Zotti, G., Granato, G. L., Maiolino, 
R., & Danese, L. 2005, MNRAS, 357, 1295 

Silva, L., Granato, G. L., Bressan, A., & Danese, 
L. 1998, ApJ, 509, 103 

Soifer, B. T., Neugebauer, G., & Houck, J. R. 
1987, ARA&A, 25, 187 

Valiante, E., Lutz, D., Sturm, E., Genzel, R., & 
Chapin, E. L. 2009, ApJ, 701, 1814 


Vignali, C., Pozzi, F., Fritz, J., et al. 2009, MN¬ 
RAS, 395, 2189 (V09) 

Yates, R. M., Henriques, B., Thomas, P. A., et al. 
2013, MNRAS, 435, 3500 

Zhao, D. H., Mo, H. J., Jing, Y. P., & Borner, G. 
2003, MNRAS, 339, 12 


This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 


26 



